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Learning Aims: You will learn to 



State the formulae for average kinetic energy of atoms, total 
energy and heat capacity of the ideal gas. 

State the quantisation condition for a particle in a cubic box. 
Explain and derive density of states for an ideal gas. 

Explain why the ideal gas gives the wrong heat capacity for 
electrons in a metal. 

State the formula for Fermi-Dirac distribution. Derive the 
formula for Fermi energy. 

Derive the interval of energy from which electrons are excited. 
Derive the distribution of the excited electrons. 

Derive the formula for heat capacity of electrons in a metal. 
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A particle in a 3-D box 



In an ideal gas, we assume that there is little interaction 
between atoms. We shall learn how to calculate the total 
energy and other properties, and apply this to electrons, 
phonons, photons, liquid helium, etc. 

From the kinetic theory of gases, we know that the average 
energy of each atom in an ideal gas is 

1 2 3, 

-mv z = -k R T. 

2 2 B 

This means that the total energy is 

U = -Nk R T 
2 a 



and the heat capacity is 



C = ™ = 3 -Nk B . 
dT 2 
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Particle in a box 



Consider a particle in a 3-D box. One solution to the 
Schrodinger's equation is a wavefunction of the form 

4>(x) = Sin(fc^a;). 





The boundary condition that the wavefunction is zero at the 
edge of the box gives this quantisation condition in each 
direction: 

rise — ■ 
a 

where n x is a positive integer. Likewise for the y and the z. 
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Density of states. 



To find the total energy of N particles, we need to add up all 
their energies: 

U = Y, n i e i- 

where e$ is the energy of level i and is the number of 
particles at that level. 

At room temperature, many of these levels could be occupied. 
This is because for an ideal gas, spacing between energy levels 
is much smaller than the average energy of a particle. 

As a result, a plot of n.j against may look like this. 

n, 1 1 II I II III II I II in 
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Ideal gas. 



The other part of the solution to Schrodinger's equation is the 
energy of the wavefunction: 

E = 



2m 



where m is the mass of the particle and k is the wavevector 
given by 

We assume that in an ideal gas, each particle (atom, electron, 
...) has the above quantisation condition and energy. 

Then we can calculate the total energy of all particles, as well 
as many properties of the gas. 



Density of states. 

This suggests that we may approximate the graph to a curve, 
and the sum to an integral. 




Within a small interval ds, the energy e is nearly constant. If we 
know the number of particles in this interval, we can multiply 
by e to find the total. 



Because of the quantisation condition, there can only be a 
certain number of states in de. For each state, there is a 
certain probability that it would be occupied. Lets start by find 
the number of states. 
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Counting states 



The wavevectors 

_n x n 
a 

are quantised at uniform intervals in all three directions. 

we can imagine a k space in which the the x coordinate is k x , 
and so on. If we use a point to represent each state, we would 
get a lattice like this. 




Because k x , k y and k z are all positive, only the points within 
1/8 of the sphere of radius k, where k is the wavevector for e. 

Since the spacing between points is n/a, the volume associated 
with each point is (n/a) 3 . Therefore, the total number of 
states is one eighth of the sphere volume divided by (ir/a) 3 : 



10 



Counting states 



First, we find the number of all states below energy e. This 
energy is related to the wavevector k by 

h 2 k 2 



2m 



k is in turn related to the components by 

k — m 

This describes a sphere in k space. 




DOS 

We find that the total number of states below k is 

, s Vk 3 
G( » = *T^> 

where V is the volume a 3 of the box. 

Suppose k is increased by dk. Then the number of states 
increases by dG. So the density of states is 

dG(k) _ Vk 2 
dk ~ 2?r : 

We have earlier defined the DOS in terms of energy, e. 



9(k) - 



To convert the variable from k to e, we must use the method 
for probability density function: 

gs(e)de = g k (k)dk. 

The subscripts are added here to emphasise that g(e) and g(k) 
are different functions. 
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Changing variables 



Rearranging, 

de 

We can now substitute the relation 

h 2 k 2 

e = 



2m 

to find gte(e). The answer is 

. . 4-mirV ,„ .1/9 
g(e) = - fcr (2me) 1 / 2 

where we have dropped the subscript again, as is the normal 
practice. For comparison: 

9(k) = tt--. 

If you substitute e for k in g(fc), you DO NOT get g{e). g(k) 
and g{e) are DIFFERENT functions! 
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Probability that a state is occupied. 




The number of particles in an interval de is equal to the 
number of states in that interval times the probability that 
each state is occupied. From statistical mechanics, this 
probability is given by the Boltzmann distribution: 

/(e) = Aexp(-e/k B T). 
where A is an unknown constant. So the number of states is 

n(e)de = Ag(e) exp(-e/k B T)de, 

where n(e) is defined as the number of particles per unit 
energy, or the number density. 
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Changing variables 



It is also useful to have the formula for total number states 



in terms of energy. 



Since no e or fc intervals are involved, we can simply substitute 

h 2 k 2 

e 



2m 

The answer is 

G{e) = 3 ^(2me) 3 / 2 . 
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Energy distribution. 
We can find A by summing to give 

TOO 

N = A g(e) exp(-£i/k B T)de, 
JO 

where N is the total number of particles. The integral 

roc 

Z SP = J 5(e) exp(-e/k B T)ds. 

is called the partition function. So we get A = N/Zgp- 
Substituting the formula for g{e), we find 

A = N( h 2 \ 3/2 
~ V \2nmk B TJ 
Substituting this into the formula for n(e), we finally get the 
full expression: 

^) = ^372^ 1/2 -P(-/W. 
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Total energy in ideal gas 



We can now find the total energy in an ideal gas. Replace the 
sum 

U = Yu n i £ i 

by the integral 

U = n{e)ede. 
JO 

Substitute the formula for number density: 

roo 2irN 1 /o 

U= ( £ ) 1 /2 e xp(^ £ /fc J? r) £ d £ . 
Jo {-Kk B Tyi 2 

Integrating this gives the formula for the total energy that we 
know from kinetic theory of gases: 

U = -Nk R T 
2 B 

This is evidence that the Boltzmann distribution for the ideal 
gas is correct. 
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Electrons in a metal 



The Drude model was propose in 1900 by Paul Drude to 
explain the transport properties of electrons in a metal. 




This treats the electrons like an ideal gas. It explains very well 
the DC and AC conductivity in metals, the Hall effect, and 
thermal conductivity. 

However, the ideal gas formula 

C = -Nk B 
2 B 

greatly overestimates the heat capacities of electrons in a 
metal. 
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Fermi-Dirac statistics 
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Why did the electron gas model fail? 



The Boltzmann distribution for the ideal gas has an 
assumption: that it is unlikely for two atoms to occupy the 
same energy level. 



This is valid for an ideal gas at room temperature, because the 
energy levels are very close together, compared to the average 
energy of the atoms. 

Unfortunately, this is no longer true for electrons in a metal. 
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The exclusion principle 



Because of the electrons are much lighter, the energy spacings 
are also much wider. So it is quite likely that two electron 
would occupy the same state at lower energy. 

However, electrons are not allowed to occupy the same energy 
states. So they have to be stacked up from bottom to top. 



!Z 
LU 



When heated, most of the electrons are stuck - there is no 
space above to move up in energy ! 

Only those near the very top can. So the heat capacities are 
much smaller than that of an ideal gas. 
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Occupying the energy states 



For convenience, we define the following function: 



g{e) exp((- M + e)/fc B T) + 1 
This is the probability that a state is occupied. 

/(e) is also called the occupation number. We need to see 
what this looks like and how it changes with temperature. 

Lets start with the simplest case: T = 0 K. If we allow T to 
approach zero, we find: 

/(e) = 1 for e < fi 
/(e) = 0 for e > n 

This means that all states with energy below fi are fully 
occupied. All states with energy above are empty. 
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Fermi-Dirac statistics 



For these reasons, we cannot use the Boltzmann distribution 
for the ideal gas. 



Instead, we need to use the Fermi-Dirac distribution: 

niejae = — — . 

exp((e - n)/k B T) + 1 

where n is an unknown constant, called chemical potential. 

The total number of electrons is fixed and is given by 
integrating the normalisation condition: 

N _ f°° ff(e)tfe 

Jo exp((e - n)/k B T) + 1 

In principle, we could solve for fi. In practice, this is usually 
difficult. 
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Occupying the energy states 



The graph for /(e) looks like this 
f(s) 



The shape shows that the energy levels are occupied below a 
certain energy, and unoccupied above that. 



This feature is characteristic of the Fermi-Dirac distribution 
that we are studying. 
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The Fermi energy 



At T = OK, The highest energy in the stack of electrons is fx. 
This energy is also called the Fermi energy, Ep. 



We can find the Fermi energy E F by integrating the 
normalisation condition: 



N = 2 x 



roo 

/ n(e)de 
JO 



The factor of 2 must be added because each energy state can 
be occupied by 2 electrons - spin up and spin down. 

We shall solve this for the Fermi energy Ep. 
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The Fermi energy 

We have previously derived the density of states: 

4rmrV , 



9(e) 



h 3 



-(2me) 1 / 2 



In the topic on ideal gas, this is obtained by counting the 
number of energy states of a particle in a 3-D box. In this topic 
on electrons, we have used the same particle in a box model. 



So the same formula for the density of states can be used for 
both the ideal gas and the electrons. We can therefore 
substitute the formula into the normalisation integral 

rE F 

N = 2 x / g(e)de 
and solve for the Fermi energy. The result is 

h 2 {3n 2 N\ 2/3 



E f 



2m 



V 
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The Fermi energy 



In terms of the occupation number, 

n(e) = g(e)f(e) 
So the normalisation condition can be written as: 

N = 2 x J n(e)de = 2 x J g(e)f(e)de 

We know that at OK, /(e) = 0 for e > fi. So So the integration 
would stop at e = fi: 

N = 2 x J^ F g(e)f(e)de 

since (j, = E F at OK. 

We also know that at OK, /(e) = 1 for e < fi. So 

rE F 

N = 2 x / g(e)de 
JO 
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Electronic heat capacity 
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Heat capacity 



We want to use the ideas and formulae that we have developed 
to calculate the heat capacity of electrons. 

To see how to do this, recall that the Drude model would 
predict for electrons a heat capacity that is the same as that of 
an ideal gas: 

3. 



C 



-Nk 



Li- 



lt is known from experiments that the actual heat capacity of 
the electrons is much smaller. 



Heat capacity 



Assume that the thermal energy is indeed much smaller than 
the Fermi energy. We shall derive an expression for this thermal 
energy, and then calculate it at room temperature to see if the 
assumption can be justified. 

At a temperature above 0 K, the occupation number 

f(e) _ 1 

exp((- M + e)/fc B r) + 1 
no longer has a sharp step at the Fermi energy. 



f(e) 



If k B T is much smaller than /x, the graph would remain close to 
the step. 
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Heat capacity 



We can understand this if we think that the electrons that are 
stacked up to the Fermi energy do not have enough energy to 
jump out of the stack. 



- 

0 
c 

LU 



This would only be true if the thermal energy is much smaller 
than the Fermi energy. 

Lets find out if this is true. 
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Heat capacity 

The "smoothened" slope of the graph tells us that electrons 
just below the Fermi energy = E F ) is excited above it. 



f(B) 



We can estimate the gain in thermal energy of the excited 
electrons from the width 8 of this slope. 



Lets take a closer look at the occupation number 

f(e) 



exp[( £ - M )/fc B T] + 1 



31 



Heat capacity 



For energy e lower than /i by a few times of k B T, The 
exponential function exp[(e — ^)/k B T] would quickly become 
small. 

The occupation number is then approximately 

exp[(e - n)/k B T] + 1 
where we have used the binomial expansion and kept only the 
first order term. 




1-exp[(s-n)/k B T] 



Heat capacity 



When temperature increases above 0 K, the step in the 
Fermi-Dirac distribution becomes smoother as electrons just 
below the Fermi energy are excited above it. 




Notice that the "tail end" of the distribution - to the right - 
looks exponential. Maybe, these electrons can behave like the 
ideal gas and approximately obey the Boltzmann distribution. 

Lets find out. 
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Heat capacity 



This means that the part of the graph to the left of fj, tends to 
the line /(e) = 1 exponentially. 




It reaches within /(e) = 1 by a fraction of 1/e, over an energy 
range of k B T. 

This means that it is the electrons within this energy range 
that is excited. So the thermal energy of the excited electrons 
is of the order of k B T. 
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Heat capacity 

For energy e higher than n by a few times of k B T, The 
exponential function exp[(e — fi)/k B T] would quickly become 
large. 



Then 



exp[(e - fi)/k B T] + 1 exp[(e - n)/k B T] 



which is just the exponential function with negative argument 
/(e) =exp[-(e-/x)/fe B T]. 




exp[-(s-n)/k B T] 
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Heat capacity 



This means that the part of the graph to the right of fx falls off 
exponentially. 

The excited electrons follow a Boltzmann distribution! 



f(e) 




Note also that the graph falls by a fraction of 1/e over an 
energy range of k B T. 
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Heat capacity 



and the number density is 

N 

— = N A + 23.71 

V A 

where N A is Avogadro's number. This gives an answer of 
2.54 x 10 28 m" 3 . 



Using the Fermi energy formula, 

h 2 {3n 2 N\ 2/3 



2m \ V 

where m is the mass of the electron, we can find that the 
Fermi energy is 3.16 eV. 
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Heat capacity 



At this point, we should justify our assumption that k B T at 
room temerature is much smaller than the Fermi energy 
H = E F , which is defined at 0 K. 

We shall take sodium metal as an example, and calculate k B T 
and E F for this metal. 

In sodium, each atom has one valence electron. This electron 
is mobile and forms the electron gas that we are talking about. 
In order to calculate the Fermi energy, we need the number 
density N/V. We can calculate this from the following data: 

density = 0.97 g cm -3 

relative atomic mass = 23.0 

So the volume for one mole of atoms is 

23^0.97 = 23.71 cm 3 , 
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Heat capacity 



In contrast, at room temperature 298 K, we can calculate that 



This is about 120 times smaller than the Fermi energy. 

We can repeat this for other typical metals, and we would get 
similar answers. 

This justifies our assumption that at room temperature k B T is 
much smaller than the Fermi energy. 

We are now a step closer to estimating the electronic heat 
capacity. Next, we need to understand the behaviour of the 
excited electrons. 
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Heat capacity 



It is clear now that at room temperature, only electrons close 
to the Fermi energy are excited. Most of the electrons are 
below the Fermi energy and are not excited at all. Since they 
cannot be excited, these electrons would not contribute to the 
heat capacity . 




So it is mainly electrons close to Fermi energy that would 
contribute to the heat capacity. We can use these to estimate 
the heat capacity and ignore the rest. 
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Heat capacity 




We know that the number of particles in a given energy 
interval is 

n(e)d£ = 2g(e)f(e)de 

where the factor of 2 again comes from the spin states of the 
electrons. 

At 0 K, the energy states below Ep are fully occupied, i.e. 
/(e) = 1. At temperature T, the number would be 

N X w 2g{E F )k B T 
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Heat capacity 



As we have seen, excited electrons behave like the ideal gas. 
We know that the energy of ideal gas is 

U X = \Nik B T 
where N\ is the number of particles in the ideal gas. 

In the case of the electrons in metal, N% refers to the number 
of electrons above the Fermi energy, and not the total number 
of electrons. We can estimate this number as follows. 
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Heat capacity 

From the definitions of the density of state and the Fermi 
energy, it can be shown that 

, . 3iV 

g(E F ) 



where N is the total number of electrons. This gives 

f 3N\ 3 k R T 

Ni « 2 k B T = -N-^— 

\AE F ) M 2 E F 

Substituting into the energy for ideal gas 

Ui = \NxkgT 

we get the energy for the excited electrons 

U X = 3 - (*N*£) k B T = 9 -Nk B k * T2 



2\2 E F J 4 E F 

Differentiating with respect to T, we get the electronic heat 
capacity 

C = lNk B ksT 



2 " E 



F 
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Heat capacity of a metal 



We have obtained the electronic heat capacity 

9 k R T 

°=2 NkB it 

More detailed calculations show that the factor of 9/2 should 
really be tt 2 /2: 

7T 2 k B T 
C = —Nk B 1 



2 E F 

For 1 mole of electrons, N = N^, and the equation can be 
written as 

2 Tp 

where T F = E F /k B is called the Fermi temperature. 
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Heat capacity of a metal 



In 1955, William Corak and his fellow co-workers measured the 
heat capacities of copper, gold and silver from IK to 5K in 
their laboratory in Pittsburgh. 

They plotted cy/T against T 2 , and got the straight lines. The 
picture here is a sketch of their results. This shows that the 
predictions of the Fermi-Dirac statistics are correct. 



Heat capacity of silver 




5 10 



T 2 (K 2 ) 
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Heat capacity of a metal 



Notice that the heat capacity is directly proportional to T. 

n 2 T 
C = —Nk B — 

This is often written in the form 

C = -yT. 

There is another contribution to heat capacity of metal. This 
comes from the vibrations of the atoms, and it is proportional 
to T 3 . 

We could imagine writing the total heat capacity in the form: 

c v = 7 T + AT 3 . 

We can measure this heat capacity to check if the formula is 
correct. Suppose that we have obtained a table of values for T 
and cy. To check if the formula is correct, we can rewrite it in 
this form: 

cy/T = 7 + AT 2 
45 



Worked Examples 
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Example 1 

A bottle contains 2 moles of neon gas at 20 oC. 

(i) Write down the formula for average kinetic energy of 
particles in an ideal gas. What is the average kinetic energy of 
the neon atoms? 

(ii) Write down the formula for total number of states. What is 
the total number of states below this energy? Assume that 
volume of the gas is 23.8 dm -3 , and relative atomic mass of 
neon is 20. 

(Mi) Write down the formula for density of states. What is the 
density of states at this energy? 
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Mass of neon atom = 20m u = 20 xl.66 x 10~ 27 kg. 



Substituting, we find that total number of states below e is 
G(e) = 2.765 x 10 30 . 



(iii) Formula for density of states is 

g(e) = — p— (2me) 1 /2 
Substituting, we find that density of states at e is 
g(e) = 6.837 x 10 50 J" 1 . 
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Solutions 



(i) Formula for average kinetic energy of particles in an ideal 
gas is 

e = -k B T. 
where T is the temperature. 

Temperature T = 273 + 20 = 293 K. Average kinetic energy 
of the neon atoms is 

e = - x 1.38 x 10" 23 x 293 = 6.065 x 10" 21 J. 
2 

(ii) Formula for total number of state: 

OW = ^( 2 -) 3/2 - 
Volume V = 23.8 dm" 3 = 23.8xl0~ 3 m~ 3 . 
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Example 2 

A piece of copper contains one mole of atoms. Each atom 
gives one conduction electron. 

Suppose that these electrons form an ideal gas. 

(i) Find the total energy of these electrons at 300 K. 

(ii) Find the average kinetic energies of the electrons. 

(iii) Find the heat capacity of the electrons. 
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Solutions 

(i) Formula for total energy of one mole of ideal gas is 

3 

U = -RT. 
2 

where T is temperature. Given T = 300 K. So 

U = - x 8.31 x 300 = 3739.5 J. 
2 

(ii) Formula for average kinetic energy of particles in an ideal 
gas is 

e = |fc B r. 
where T is the temperature. 

Average kinetic energy of the electrons is 

e = - x 1.38 x 10" 23 x 300 = 6.210 x 10" 21 J. 

2 
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Example 3 

In fact, only the excited electrons in copper behave like an ideal 
gas. 

(i) Write down the correct formula for the heat capacity of the 
electrons in terms of the Fermi energy. 

(ii) The volume for one mole of copper is 7.11 cm 3 . Each 
copper atom gives one conduction electron. Find the Fermi 
energy. 

(iii) Find the correct heat capacity of the copper at 300 K. 
What percentage is this compared to the ideal gas result in 
Example 2? 
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(iii) Heat capacity of the electrons is 

C = -R = - x 8.31 = 12.47 J/K. 
2 2 
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Solutions 

(i) The correct formula for heat capacity of the electrons is 

(ii) Fermi energy is 

E F 



h 2 {3n 2 N^ 2/3 



2m \ V 



m = m e 
N = N A 

V = 7.11 cm 3 = 7.11 m 6 . 



Substituting, we find 



E F = 1.117 x 10" 18 J. 
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(iii) Substituting the above values, the correct heat capacity of 
the copper is 

C = -NkJ^- = 0.1386 J/K. 
2 a E F 

Compared with the ideal gas calculation, this is 
0.1386/12.47 x 100 = 1.112%. 
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Solutions 



(i) Formula for total number of states: 



(ii) Two electrons can occupy each energy state. So total 
number of electrons states below the energy e is 2G(e). 

At 0 K, all electrons are below E F , so total number of electrons 
is: 

N = 2G(E F ). 
Substituting the total number of states formula: 

N = 2 x fj(2m^) 3 / 2 
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Example 4 



(i) Write down the formula for total number of of states, G(e), 
below energy e for particles in an ideal gas. 

(ii) Derive the Fermi energy formula using G(s). 

(iii) Derive g(E F ) in terms of E F and number of particles N 
only (without V, m or h). 
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Rearranging: 



r? {3n 2 N\ 2/3 



2m \ V 



(iii) Notice that the total number of states may be written in 
this form: 

G(e) = Ke 3 ' 2 . 
where K is a factor that does not depend on e. 

Taking log on both sides and differentiating with respect to e, 
we find the density of states 

g(e) = 3 -Ks^. 
Dividing this by the previous equation, we eliminate K: 

9(e) = 3 

G(e) 2e' 
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At Fermi level, e is E F : 

9(E F ) _ 3 
G(E F *) 2i?p 

Since G(E F ) is the total number of states below Fermi level, so 
the total number of electrons is 

N = 2G(E F ). 

Substituting this into the previous equation and rearranging, we 
get 

^ 3iV 
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Learning Outcome: You should be able to 



State the formulae for average kinetic energy of atoms in an 
ideal gas, the total energy and heat capacity of the ideal gas. 

State the quantisation condition for a particle in a cubic box. 
Explain and derive density of states for an ideal gas. 

Explain why the ideal gas gives the wrong heat capacity for 
electrons in a metal. 

State the formula for Fermi-Dirac distribution. Derive the 
formula for Fermi energy. 

Derive the interval of energy from which electrons are excited. 
Derive the distribution of the excited electrons. 

Derive the formula for heat capacity of electrons in a metal. 
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